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Abstract 

In the present article we consider some general class of divergent 
diagrams in Keldysh diagram technique. These divergences arise for 
non-equilibrium matter and closely related to the divergences in the 
kinetic theory of gases. We suggest a scheme of renormalization of such 
divergences and illustrate it on some examples. In the other papers of 
these series we develop the general theory of renormalization of non- 
equilibrium diagram technique. The fact that thermal divergences in 
non-equilibrium diagram technique can be renormalized leads to the 
following consequence: to prove that the system tends to the thermal 
equilibrium one should to take into account the behavior of the system 
on its boundary. In this paper we illustrate this fact on Bogoliubov 
derivation of kinetic equations. 
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1 Introduction 



There are exist divergences in virial decomposition of kinetic equations. This 
fact was observed by Cohen and Dorfman [1]. It is possible to summize some 
set of diagrams to obtain a finite result. It was done by Kawasaki and 
Oppenheim [2]. 

Our main goal in this series of papers is to analyze such divergences. 
In the present paper we study divergences in Keldysh diagram technique 
which arise if the state of the matter is non-equilibrium. It is more or less 
obvious that these divergences are the same as the divergences in the kinetic 
equations. 

In the present series of papers we develop the general theory of such 
divergences. As result for a wide class of Bose systems in the sense of formal 
power series on coupling constant we find non-Gibbs state (•) such that the 
correlators 

(V ± (t,x 1 )...V ± (t,x n )) (1) 

are translation invariant, do not depend on t and satisfy to the weak cluster 
property. Here ^Z 1 * 1 are the fields operators and the weak cluster property 
means the following 

lim / (^ ± (t,xi + 5 1 e 1 a)...^f ± (t,x n + 5 n e 1 a))f(x 1 , ...,x n )d i Xi...d i x n 

|o|->oo J 

= J (* ± (t,x il )...* ± (t,x i J)(* ± (t,x ik )...* ± (t,x in )) 

xf(x 1 ,...,x n )d 3 x 1 ...d 3 x n , (2) 

there Si G {1,0}, i = l,2...n and 

ii < i 2 < ... < ik, 
ik+i < ik+2 < ■■■ < in, 
{h,i 2 , .-.,4} = {i = l,2...n\5i = 0} ^ 0, 
{4+1,^+2, i n } = {i = l, 2...n\5i = 1} ^ 0. (3) 

f(xi, ...,x n ) is a test function, e± is a unite vector parallel to the x-axis. In 
the present paper we consider only a some wide class of divergent diagrams. 

Let us prove that the existence of such states implies non-ergodic property 
of the system. We consider the problem only on classical level. Suppose that 
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our system is ergodic, i.e. there no first integrals for the system except energy. 
Then, the distribution function depends only on energy. We can represent 
the distribution function f(E) as follows: 

f(E)=Y,c a 8(E-E a ), (4) 

where the sum can be continuous (integral). Let 1 be some finite subsystem 
of our system. Let 2 be a subsystem obtained from 1 by translation on the 
vector L parallel to the x-axis of enough large length. Let 12 be a union of 
the subsystems 1 and 2. Let p±, p 2 and p i2 be distribution functions for the 
subsystems 1,2 and 12 respectively. Let r\, T 2 and be points of the phase 
spaces for the subsystem 1, 2 and 12 respectively. By the same method as 
the method used for the derivation of the Gibbs distribution we find: 

P\2 = 2_^c a — — - (5) 

Z/q, Zj a 

in the obvious notation. But the weak cluster property implies that 

P12 = PiP2- (6) 

Therefore all the coefficients c a are equal to zero except one. We find that 

f(E) = cS(E - E ) (7) 

for some constants c and E Q . So each finite subsystem of our system can be 
described by Gibbs formula and we obtain a contradiction. 

Non-ergodic property means that there no thermalization in infinite Bose- 
gas system. 

This fact implies that to prove that the system tends to thermal equi- 
librium we should to take into account the behavior of the system on its 
boundary. Indeed if the system has no boundary the system is infinite. 

To illustrate this fact we will study Bogoliubov derivation of kinetic equa- 
tions. When one derives BBGKI-chain one neglects by some boundary terms. 
If one take into account this boundary terms and use Bogoliubov method of 
derivation of kinetic equations one find that these boundary terms compen- 
sate the scattering integral. 

The paper composed as follows. In section 2 we introduce the notions of 
the algebra of canonical commutative relations. In section 3 we describe our 
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model. In section 4 we describe non-equilibrium (Keldysh) diagram tech- 
nique. In section 5 we discuss the divergences in our model and the method 
of its renormalizations. In section 6 we give a proof that there exists diver- 
gences in our theory. In section 7 we describe regularization which will be 
used. In section 8 we discuss some simple relation on Green functions. In 
section 9 we make renormalization procedure in one-chain approximation. In 
section 10 we begin to renormalize diagrams in two-chain approximation. We 
calculate divergent parts of all diagrams in this approximation. In section 
11 we discuss subdivergences i.e. we calculate the contributions which comes 
from counterterms for the one-chain diagram. In section 12 we show that the 
divergent part of all diagrams which is proportional to -^8{ui — u(p)) can be 
subtracted by counterterms. In section 13 we show that the divergent part 
of all diagram which is proportional to \5'{uj — co(p)) can be subtracted by 
counterterms. In section 14 we study Bogoliubov derivation of kinetic equa- 
tions and show that the scattering integral is compensated by some boundary 
terms which are usually neglected. Section 15 is a conclusion. 



2 The algebra of canonical commutative re- 
lations 



The algebra of canonical commutative relations is a * algebra with a unite 
generated by generators 



where / belongs to the Schwartz space of test functions >S(1R 3 ). The 
generators a(f), a + (f) satisfies the following relations: 
a(f) is an antilinear functional on /, 
a + (f) is a linear functional on / and 



(8) 



Hf),a + (g)] = (f,g) 
[a(f),a(g)] = [a + (f),a + (g)]=0. 



(9) 



Here (/, g) is a standard scalar product in L 2 (R 3 ) 




(10) 



4 



The involution * is denned by the following rule 

W)y = a + (f). (11) 
The "field operators" ^(x) are defined by the following formulas 

9(x) = -i-r / e ikx a(k)d 3 k, 

(27TJ2 J 

= — ^ [ e~ lkx a + {k)d 3 k. (12) 

(27r)a y 

Here we have used the following formal notation 



a+ 



(/) = f a(k)f(k)d 3 k. (13) 



Definition. Let us consider a Gauss state po on the algebra of canonical 
commutative relations defined by its two-point correlator as follows 

Po (a(k)a(k')) = Po (a + (k)a + (k>)) = 0, 

p (a + (k)a(k')) = n(k)5(k - k'). (14) 



If 



e -/M*0 

n W = i- e -Mfc) ' (15) 



the state p is called a Plank state. Here u;(fc) — \ — p, p < 0. 

3 The model 

Our model is described by the following Hamiltonian 

H = H + XV, (16) 

where 

H = J uj(k)a + {k)a(k)d 3 k, 

u(k) = ^-p, p<0, (17) 



V = \ I ^ + (x)^ + (x')V(x - x')^{x')^(x)d 3 xd 3 x', 
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AeK, (18) 



V is an interaction, V(x) £ S(M. 3 ). Let us rewrite the interaction in the 
Fourier representation 

V = 1 J (^h^hd^k'^k^Vih + k 2 ) x 
5(fci + fc 2 - ifci - k'^a+ik^a+ik'^aik^aih). (19) 
Here by definition 

V(k) = J e ikx V(x)d 3 x. (20) 

4 Nonequilibrium diagram technique 

Let us introduce the Green functions for the system 

p(T(*%(t 1 ,x 1 ),...,*%(t n ,x n ))). (21) 

Here symbol H near ^ r± means that are Heizenberg operators. 

In nonequilibrium diagram technique we admit the following representa- 
tion for the Green functions 

p(T(*%(t 1 ,x 1 ),...,*%(t n ,x n ))) = 
p(5- 1 T(*±(t 1 , xx), *±(i n , x n )5)). (22) 

Here the symbol near SE^ means that ^ are operators in the Dirac repre- 
sentation (representation of interaction). The S'-matrix has the form 

+oo 

S = Texp(-i J V(t)dt), (23) 

— oo 

and 

+oo 

S- 1 = Texp(z J V(t)dt). (24) 



Here T is a symbol of antichronological ordering. 

Let us recall the basics elements of Nonequilibrium diagram technique. 
The vertices coming from T-exponent are marked by symbol — . The vertices 
coming from T-exponent are marked by symbol +. There exists four tips of 
propagators 

G£"(ti - t 2 , Xl -x 2 ) = p (^(t 1 ,x 1 )y + (t 2 ,x 2 )), 

G + (t! - t 2 , Xl ~ X 2 ) = Po(^ + (t 2 , X 2 )V(t 1 ,X 1 )), 

G Q -{t l -t 2 ,x 1 -x 2 ) = Po (T(y(t 1 ,x 1 )y + (t 2 ,x 2 ))), 
- t 2 , x x - x 2 ) = po(TW*i, xi)^ + (h, x 2 ))). 

(25) 

Let us write the table of propagators 

G+-(t,x) = J ^(27r)%-o;^))(l + ^))e-^), 

G +(t,x) = J ^(2tt)*( W - u{k))n{k)e-^- k *\ 

C -(t r ) -i [ d ' k r ri(k) -i(ut-kx) 

G ++( t x) -i f d ' k r < k ) 1 + n(k) , M _ M { . 

5 Divergences 

A typical example of divergent diagram is pictured on figl. 
fig. 1 



The ovals 

represent the sum of one-particle irreducible diagrams. These diagrams are 
called chain diagrams. Let us suppose that all divergences of self-energy 
parts (ovals) are subtracted. The divergences arises from the fact that 
singular supports of propagators are coincide. At first we consider diagrams 
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with one self-energy insertion (one-chain diagram). These diagrams are 
pictured on fig. 2. 

fig. 2 



+V J + 



These diagrams are analogues to 

one-loop diagrams in quantum field theory. 

There exist two possible types of counterterms. The first one is a coun- 
terterms of mass renormalization. Mass renormalization is equivalent to the 
following replacement 

XV -> XV + M, (27) 

where M has the form 

M = J m(k)a + (k)a(k)<Pk, (28) 

m(k) is a real-valued function of k. 

The second type of counterterms are counterterms of asymptotical state. 
Asymptotical state renormalization is equal to the following replacement 

1 - f h(t)dt 
Po(-) - ^Po(e - (•)), (29) 

where 

h = J h(k)a + (k)a(k)d 3 k, 
h{k) is a real-valued function and 

- / h(t)dt 

Z = Po (e -~ ). (30) 

We will proof below that the counterterms of asymptotical state are enough 
for the renormalization of all one- and two-chain diagrams. 
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6 Proof of the existence of divergences in the 
theory 

Suppose that there no divergences in Keldysh diagram technique if n(k) ^ 
— ^ for any positive a, (3. Therefore the Green function 

e «T-+/5 + l 

p(S- 1 (Sy + (t 1 ,x 1 )y(t 2 ,x 2 ))) (31) 
is translation invariant. So the density matrix 

Pt(xi,x 2 ) := p(S- 1 (SV + (t,x 1 )*(t,x 2 ))) (32) 
is an integral of motion. Let 

p t (k) = J d 3 xp t (0,x)e ikx . (33) 

In zero order of perturbation theory p(k) = n{k). But if there no divergences 
in Keldysh diagram technique it is possible (see [3]) to derive the following 
kinetic equation for p(k) 



J 



dptjk) 
dt 

w{p,Pl\P2,Ps){(l +p(p))(l + P(P1))P(P2)P(P3) 

-p(p)p( Pl )(l + p(p 2 ))(l + p{ps))}. (34) 
The right hand side of this equation is equal to zero only if 

pw = -^A — ( 35 ) 

for some a, (3 (a > 0, (3 > 0). But n(k) = p{k) in zero order of perturbation 
theory so n(k) has a Bose-Einstein form. This contradiction proves our 
statement. 

7 Regular izat ion 

Let us now introduce regularization. Note that 

' ' (36) 



x + is x 2 + e 7r x 2 + e 2 
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Therefore we use the following regularization 



5{uj — oo(k)) — > — 
1 



n 



n (to - uj(k)) 2 + e 2 ) 
uj — u(k) 
oj-u(k) ) ~* {uj - uj(k)) 2 + e 2 ) ' 



--: 5 £ (uj - Lo(k)), 
1 



uj — u(k) 



8 Some simple relation on the Green 
tions 

Lemma 1. The following equalities hold 

G (ti - t 2 , xi - x 2 y = G ++ (t 2 -ti,x 2 - Xi), 
G + ~(t 1 - t 2 , x\ - x 2 )* = G + ~(t 2 - ti, x 2 - xi). 

Proof. We have 

G (ti - t 2 , x x - x 2 y = p (T(*(t 1} xi)* + (t 2 , x 2 ))Y 
= p (f(y + (t u x 1 )V(t 2 ,x 2 ))) = G ++ (t 2 - t u x 2 -x x ). 

So the equality [38] is proved. We have 

G + -(h - t 2 , x x - x 2 y = p (*(ti,xi)* + (t 2 , x 2 )Y 
= PD(*(*2,z 2 )tf + (*i,xi)) = G + -(t 2 -t u x 2 - x x ). 

So the equality EH is proved. 
The Lemma is proved. 
It is easy to prove the following 
Lemma 2. The following equality holds 

G+-(t, x) = 9(t)G—(t, x) + 9(-t)G ++ (t, x), 
G- + (t,x) = 9{t)G ++ {t,x) + 9(-t)G—(t,x) 

Let us introduce the following matrix 



G 



C++ G H 
G~+ G~ 
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Let us introduce the similar matrix for the self-energy operator 

E = 



E+ H 



E~ 
E" 



Dyson equations in Fourier representation has the form 

G — Go Gq,EG. 
We have from these equations that 



G 1 



C 



-i 



or in the matrix form 



detG 



G 
— Gr 



_ C(3 

G n 



detG 



G~ 
-G- 



— G H 
G+ 



It follows from Lemma 1 that 



G ++ (u,p) = G-(uj,py 
G + -(cu 1 p) = G+-(u, P y i 
G- + (u,p) = G- + (u, P y. 

Therefore detG , detG are real and we have the following 
Lemma 3. 



e (ti - 1 2 , xi - x 2 y 

E +_ (ti - t 2 ,x t - x 2 ) 

The following Lemma holds. 
Lemma 4. 

E++(u;,p) + E-(^) 



E +H 

:E + 



(t 2 — t\,X 2 — Xi), 
"(* 2 - t U X 2 - Xi). 



Proof. The statement of lemma follows from the Dyson equat 
and the following two obvious equalities: 

G ++ {uo,p) + G-(u,p) = G- + {uo lP ) + G + -(cu,p), 



G+ + (uj, P ) + G (uj,p) 



G + (u,p) + G+ 



[CO 



,p)- 
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9 Calculation of the propagators in one-chain 
approximation 

Lemma 5. The following limit equality holds (in the sense of distributions): 

lim$(s) - -L* e (s)) = 0, 

e^O AITS 
1 1 

hm(-6 £ (x) - -5(x)) = reg, 

lim{- 1 - -S(x)} = reg, 

Urn 5 £ (x)V £ {-)reg =, 

xeR. (53) 

Here reg means some correct distribution. 

Proof. Let f(x) be some test function with compact support. We have 
J S 2 £ (x)f(x) = ^J ^ 2 ] r£2)2 f{x)dx 

= ^ / ^4w {/(0) + xf,{0) + x2 ^ x)}dx (54) 

for some smooth bounded function ip(x). We have 

=^\Sj^w dx)m+0 ^ (55) 



But 



So 



Sw^w dx= l (56) 

/ ^(x)f(x) = ^f(0) + O(e). (57) 
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So the first equality is proved. One can prove other three equality by the 
same way. 

Therefore we see from the Lemmas 1,2, that we can only consider the 
function G (t,x). But the function G (t,x) can be represented as a sum 
of chain diagrams. At first let us consider the diagrams with one self-energy 
insertion (one-chain diagram). We have G~~ = H % J ■> where the diagrams 

i,j=± 

for H l J are presented at the fig. 2. We have the following representation for 
the divergent parts of these diagrams. 

(H~~) div (uj,p) + (H+ + ) div (cu,p) 

= 2^--{^p)n{p){l + n{p))\8{u-u{p)) 

+2nZ ++ (uj,p)n( y p)(l + n(p))^6(u - u(p)). (58) 

We see that the divergent part of these two diagrams is real (because E = 
(£++)*). 

It is obvious from previous calculations that the sum of two possible mass- 
renormalization diagrams is equal to zero. Let us consider the singular part 
of other two diagram presented at fig 3. 

fig. 3 




We have 

(H~ + ) div (uj,p) + (H+~) div (u;,p) 
= 7r(27r)(l + 2n{p)){\ + n(p))5 2 £ (uj - u;(p))£-+(u;,p) 
+7r(27r)(l + 2n{p))n{p)5 2 e {uj - uj{p))Z+-(u,p) 

= tt(1 + 2n(p))(l + n(p))-6(u - u(p))^- + (u,p) 

e 

+ .(l + 2n W )n(p)i^-.(p))E + -( W ,p) + 0(£ ). (59) 

We see that (H~ ) div {uJ,p) + (H+ + ) div (uj,p), (H- + ) div (uj,p) and 
(H+~) div (u;,p) are real. 
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We will use the dotted line for lines which connects creation-annihilation 
operators with operators coming from the vertex: f h(k)a + (k)a(k)d 3 k (see 

fig. 4 



+v J + 



fig. 4). So the divergences in G £ 

H l J can be subtracted by the following counterterm: 

i,j=± 



h{p) = Z ++ (uj,p) + E—(u,p) + 



2n(p)(l + 

x{(l + n(p))E _+ (w,p) + n(p)£+~(^,p)}. (60) 



We have by using Lemma 4 



Wri = 1 + 2nip) x 
W 2(1 + ra(p))ra(p) 



{(1 + n{p))YT + (u,p) - n(p)Z + -(u,p)}. (61) 

The left hand side of this equation can be rewritten as follows (in approxi- 
mation used in [3]) 

- dw^m st(Ph (62) 

where St(p) is a scattering integral. So h(p) ^ for non-equilibrium matter. 



10 Calculation of propagators in two-chain 
approximation 

We will calculate the divergent parts of all diagrams H l J kl , i,j,k,l = ± 
presented at fig. 5. 
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We have 

= Yl H tj + Yl H l J kl in two-chain approximation. 



i,j=± i,j,k,l=± 



Let us calculate H e ~~ (uj,p) (fig. 6). We have 

1 + n(p) n(p) 



H-~ -M = -(E-fx l 



lo — u(p) + ie oo — u(p) — 26 



) 2 {- (i+Mp)) 3 Mp)) 3 



(c<j — c<j(p) + ie) 3 (a; — u(p) — ie) 3 
3(1 + n(p))n(p) 2 3(1 + n(p)) 2 n(p) 



1 ffi*^ 

(a; — u(p) + ie)(u — u(p) — ie) 2 (uo — uo(p) + ie) 2 (uo — uo(p) — ie) 

But ( jjz^j±^ j 1S a distribution. So we have the following expression 
for the singular part of H £ . 

3in( P ) 2 (l+n(p)) 



(a; — a;(p) + ze)(a; — c<j(p) — ie) 2 



3m(A;)(l + n(fc)) 2 
+ (a;-a;(p)+ ? £) 2 (cu-^(p)-2£) + 



We have 



(c<j — cj(p) + ie) 2 (cj — uj(p) — ie) 



(a; — a;(p) + ie) 2 (a; — uj(p) + ie) 

97T7 

= 7 " . + (65) 

(a; — a;(p) + ?e) 2 

Let J(c^) be a test function /(a;) G S". We will calculate 

J ' := ^/ (ii lwi /(gii)diI (66) 
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But f(euj) = /(O) + eQf'(0) + .... Substituting this decomposition into last 
equation, we find 

+ h m Jwrim^j dn+oir >- < 67) 

We use the Cauchy theorem for calculation these two integrals. Let us close 
the integration contour in the upper half-plane. The integrand has only one 
pole at the upper half-plane at the point Q = i. Therefore 

+00 

—00 

= 2 "(?77^ = 2 " x =5 = ~T (68) 

By the same way we find 

r n 

B:= J (n + i )(w + i) dQ 

= 2 "wnr^--T < 69) 



So we have 



or 



Ie = -^\f(0)-^m + O(l), (70) 



2 

5 £ (lo - u{p)) 



uj — uj(p) + is 



1 i 



+0(1). (71) 



In result 



1 



(a; — u>(p) + ie) 2 (u — u(p) — is) 
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2*i(-^*( W " "(p)) + - w(p)) + 0(1)). 



(a; — u;(p) + 2e)(c<j — uj(p) — ie)'' 



-2m{~8^ - u(p)) - ^5'{u - u(p)) + 0(1)). (72) 



At last 



37T 

{H--~-)<Uv(w,p) = {^n(p)(l + n(p))(l + 2n(p))<y(w - w(p)) 

' n(p)(l + n(p))5'(c-c(p))}(S-(c,p)) 2 . (73) 



4e 

Let us now calculate the diagrams presented at fig. 6,7. 
fig. 6 



+V J + 



fig. 7 



+v J + 



way as previous we find 



— n(p)(l + n(p)){- 



-ra(p)(l + n(p)){-(l + 2n(p))5(w - w(p)) 



Let us now consider diagrams presented at fig. 8, 9. 



By the same 



--5'(co - cu(p))}Z-(u;,p)Z ++ (u;,p). (74) 
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fig. 9 



+ 



It is easy to 

see that 

H-+-(u,p) = (1 + y ) S- + (a;,p)(S + -(a;,p))- 1 g £ -- + "(^ ) P)- (75) 

Omitting the calculations we find 

(H- + —) div (u;,p) = (H- + —) div (uj,p) 
= 2 7 r(l + n(p))S^ + (cu,p)S"(a;,p) 

x{[h(l + n{p)f + n{pf) + n(p)(l + n(p))} ^ ~ " {p)) 
4 e z 

+J(l + 2n(pM W - W (p))} (76) 

and 

(^7 _ + _ )^(^,P) = (#£ + )dw(^,P) 

= 27rn(p)S + ~(cj,p)S (u;,p) 
x{[i((l + n{p)f + n(p) 2 ) + n(p)(l + n(p))] ^ 

+^(l + 2 n (p))%- W (p))}. (77) 
Let us now present analytical expression for other diagrams: 

(# £ ++++ W(^P) = {(l + n(p))n(p)(l + 2n(p))^(a;- W ( P )) 

+-<5'(u; - w (p))(l + n(p))n(p)}(S ++ (c,p)) 2 , (78) 
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x(l + 3n(p) + 3n(p) 2 )n(p)5(uj — u)(p)), 

{H; +++ ) dm {u,p) = Z~ + (u,p)Z+ + (u,p) 
x(l + n(p))(l + 3n(p) + 3n(p) 2 )^5(u - u{p)), 

(i/+- ++ )^(^,p) = 37rS + -(cu,p)S ++ (cu,p) ( 5(a;-a;(p)) 

(1 + n(p))n(p) 2 



x- 



£ 2 



( J ff+ + - + )^(^,p) = 37rS ++ (^,p)S- + (cu,p) ( 5(a;-a;(j9)) 

(1 + n(p)) 2 n(p) 



x- 



£ 2 



27T 

= ( T )n(p)(l + n(p))£+-(c,p)£->,p) 
xra(p)(l + n(p)){-^(l + 2n(p))6(u - u(p)) - l -5\u - cu(p))}, 

2ix 

{H; ++ -) div {u, P ) = ^~ 2 Mp) 3 + Mp)(i + n( P )) 2 

+2n{j>) 2 {\ + n{p)) + (1 + n(p)f} 
x^- + (u,p)i: + -(u,p){l + n(p) + n(p) 2 }, 

2tT 

= --(1 + n(p)) 2 S- + (c,p)S- + (u;,p) 
x{-4(l + 2n(p))5(u - W (p)) + V(w - W (p))}, 

£ z £ 
27T 

= - T (n(p)) 2 S+-(u;,p)S + -(a;,p) 
x{-|(l + 2n(p))*(w - u;(p)) + V(w - w(p))}. 
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11 Counterterm diagrams 

Let us recall that we renormalize the asymptotical state p by the following 

way 

+00 

- f h(t)dt 1 

(87) 

where 

- / fc(t)dt 

Z = p(e - ), (88) 

and 

h{t) = e ltHa he- itHo . (89) 

So we have to take into account the counterterm diagrams pictured at fig. 
10, 11. i.e. # 00±± and H ±±0 ° respectively. 

fig. 10 



OV^VO +v J± 



fig. 11 



+ 



Let us calcu- 
late H°°— . We have 

{Hf-) div = -h(cu(p),p)(2n) 2 5 2 £ (u; - u{p))n{p)(l + n(p)) 

xz£ (u;,p){ — : — -} 

uj — uo[p) + xe uj — u[p) — is 

= -i(2ir) 2 h(p,u}(p))Y,~-(u,p)n(p)(l + n{pj) x 
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-n(p)[^)6(u - u(p)) - ~6'(u - u(p))}}. (90) 

In result 

{Hf-) div {u,p) = -h(co(p),p)E-(co,p)n(p)(l + n(p))(2n) 2 

x{l|-(l + 2n(p))*( W - co(p)) - \^5\u - uj{p))}. (91) 

Let us present now the analytical expression for G 00+ ~. Omitting some 
calculation we have: 

{Hf + ~) dlv {u,p) = -h(u(p),p)(2n) 2 Z + -(p,Lj)n(p)(l + n(p)) 

x{±6(uj - w(p))^(l + 2n(p)) - 1-U'(a; - u;(p))}. (92) 

Let us present now without calculations all other counterterm diagrams: 

^>7T 

(H™- + ) div (",p) = -h(u(p),p)i:- + (u,p)-5(u-u(p)) 

xn(p)(l+n(p)) 2 , (93) 

(#° 0+ -%>,p) = --h(u(p),p)Z++(u,p) 

xn(p)(l + n(p)) 2 5(uj - u(p)), (94) 

(H;-°°) dzv (uj,p) = -h(u(p),p)E—(u},p)n(p)(l + n(p)) 
x(2n) 2 {^8(u; - w(p))^(l + 2n(p)) - ~^(a; - a,(p))}, (95) 

(tf £ - + °V(cc>,p) = -/i(a;(p),p)S- + (a;,p)n(p)(l + n(p)) 
x(2vr) 2 {l^ - w(p))^(l + 2n(p)) - 1-U'(u, - W (p))}, (96) 

(fl^W^P) = --n(p) 2 (l + n(p))E +± Ma;(p),p) 

x«5(w-w(p)). (97) 

We have presented all counterterm diagrams. 
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12 Annihilation of all strong divergences in 
the Green functions 

Let f(e) be a function of e of the form 

f(e) = c 2 \ + Cl - + 0(l), ase^o. (98) 



e 2 e 



Put by definition 



(/(e)) 2 = c 2 



2 : 



= ci 1 -. (99) 

We find that ++ ) 2 ,(#++— ) 2 , (#- ++ -) 2 ,(tf+- + -) 2 ,(# £ -+-+) 2 , 

(H^ h ) 2 are real. But the following terms (if^ ) 2 and (i7+ +++ ) 2 are 
complex-conjugated to each other. 

All not real counterterm diagrams have the form: 

(H?—(u,p)) 2 = -(27r) 2 h(u(p),p)n(p)(l + n{p)){\ + 2n{p)) 

3 1 

x- — ^—(uj,p)S(u; -u(p)), 

(tf £ 00 ++(c,p)) 2 = -3nh(uj(p),p)n(p)(l + n(p)f 

x^ + (u,p)5(uj-u;(p)), 

(Hr °^,p)) 2 = -h(u(p),p)n(p)(l + n(p))(l + 2n(p))^-(27r) 2 

xT,~~(u,p)S(uj - u(p)), 

3tt 

7 2 

xZ ++ (u,p)5(uj-uj(p)). (100) 



[H+ +00 (u;,p)) 2 = -h(u;(p),p)n(p) 2 (l + n(p)) 2 



So it is easy to see that the sum of counterterm diagram is real. 
One can see that the sum 

{Ht + ~ + ) 2 + (# £ - +++ ) 2 + {h; + -) 2 + (Hr~ + h (ioi) 

is real and 

(Ht ++ 'h + m- ++ ) 2 + (H+—h + (Hr + -h ( 102 ) 

is real too. So all the most strong divergences can be renormalized by renor- 
malization of the asymptotical state. 
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13 Annihilation of all weak divergences in the 
Green functions 

Now we try to answer the question: if the divergences which are proportional 
to 8'{<jj — u>(p)) in Green functions are vanished. 
Let us recall what 

£-+ + £+- = -S++-E— . (103) 
Here and below we omit arguments (u>,p) of all functions. We have 

(H-—)i = -(E-) 2 ^(l + nK- (104) 

Here and below we will omit an argument of 5-function. Corresponding 
counterterm is equal 



(H- c — ) l = (Z-f-n(l + n)5'. (105) 



Therefore 



(H--- ) 1 = -(E--) 2 gn(l + nK. (106) 

Here we put by definition 

Hf£±± = + H^. (107) 

We have also 

(#++++)! = (£++) 2 ^n(l + n)5', (108) 

and corresponding counterterms are equal to zero. In result 

(H; R — )x + (^ +++ )i = -| x {(£ ) 2 - (E++) 2 }n(l + „)*'. 

This quantity is real. 

Now let us find counterterms to the diagrams 

H+ + — and H~- ++ . (109) 
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By using the table of divergent parts of diagrams we find 

(fl^-) 1 = +^n(l + n)E-E++. 
By the same way we find 

(^c ++ )i = +|^(l+^)S-S ++ , 

so 

Now we will calculate the diagram H~ + ~ + . 

W + - + )i = -£(i + ^(s- + )V. 

It is easy to find that its counterterm is equal to 

(^ + )i = ^(l + 2n)(l + n)(S-+) 2 5'. 



Therefore 



(# £ - + - + )i = -£(i + -)(£- + )V. 



Now let us calculate the diagram H + ~ 



and 



In result: 



Ht c - + ~) 1 = ^{l + 2n)n{i:+-Y5> 



Now let us consider the diagram: 

— 2ir? 

(H- + -h = — (1 + n(p))(l + 2n(p))S- + S- 
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Let us find counterterms to this diagram: 



Therefore 



^ + 1i = ^- + r(i + nH' 

+(l + n)(l + 2n)—E- + E—S' 
Ae 

= -^(l + n)I!- + I!— ^ 
Ae 



= + n)(3 + 8n)*'. 

Let us consider the diagram h : 

(i/ £ — + )i = -^(1 + + 2n)S- + S-"5'. 
The counterterm corresponding to this diagram is equal to 

= |(1 + n)(l + 2n)E- + E"5'. 

In result 

W)i + " + )i = ^(1 + n) 2 E-+E"5'. 

Now let us calculate counterterms to the following diagrams: 

{Ht— )i + 

= n (i + 2n)<J'. 



, n(l + 2n) 
4e 
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Therefore 



and 



— 1 7T 



It is easy to see that: 



(h: c = o, 



So 



+ = + n)5>. 

It is easy to find that 

(i/ £ - +++ )i + (^ £ ++ - + )i = o, 

(i^+) 1 = n(l + n)|E-+E++^ 



Therefore 

{H;+ ++ + 1Z++- + )i = n(l + n)g5T+E++*'. 
Now we must to calculate the following two diagrams: 

and {H;++-\. 
Let us start with (Hf + )i. We have 
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Therefore 

2m 

Let us calculate the diagram H~ . We have 



H+~- + ) x = - — E-+E+-n(p)(l + n(p))<y'. (136) 



In result 



#; ++ - = -^{1 + n + n 2 }E- + E + -5', (137) 



(fl" e -+ + -) 1 = S- + S + -^((l + 2n)(l + 2n))<y'. (138) 



TOT" + ^" + ) = (139) 

Now we must summarize all these contribution neglecting by real parts. 
All real parts can be subtracted by counterterms of asymptotical state. We 
have 

E 

i,j,k,l=± 

= {-j(l + -)(S- + ) 2 + ^n(S+-) 2 

-^E- + E"(l + nf - ^E+-E-n 2 

+ -E ++ E + -ri(l + n) + — + n)E-+E ++ 
2e 2s 

-£e+"E-+}5'. (140) 

Let us unite in the r.h.s. of the last formula 3rd and 6th terms, and 4th and 
5th terms. Neglecting by some real terms we find: 

E ( H st l )i 

i,j,k,l=± 

= {-£(£- + ) 2 + ^n(E+-) 2 

--E-+E ++ (l + n) + -E ++ E+-n 
2e 2e 

-— E+-E-+15'. (141) 
As 
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It follows from the identity 

£++ + IT" = -E~ + - E+~ (142) 

that 

+ = -E~~ (143) 

and 

E ++ + E + ~ = -E - (144) 

Let us unite in (11411) the first term with 3rd term and second term with 4th 
term. We find 



i ,j,k,l=± 



ijkl\ 
eR H 



= _ E — - E+~)(l + n)<5' 

+-E + -(£++-E— -£-+W 

4 

(145) 

But X + ~ and £~ + are real and (£ )* = Neglecting in (1145ft by real 

terms we find 



i,j,k,l=± 

= {^S+-S-+(l + n) - — E+-E-+n 

-^E + -E-+}5' = 0. (146) 

So the image part of divergences of Green function is equal to zero. There- 
fore Divergences of two-chain diagram can be subtracted by counterterms of 
asymptotical state. 



14 Notes on Bogoliubov derivation of Boltz- 
man equations 

In this section we study the problem of boundary conditions in Bogoliubov 
derivation of kinetic equations [4]. Let us consider N particle in R 3 . Let be 
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a coordinates of particle number i, and Pi be a momenta of particle number 
i, i = 1, N. Suppose that particles interacts by means the pair potential 
V(qi-qj). We suppose that V belongs to the Schwartz space. Let X{ = (j>i, 
be a point in the phase space T. Let f(xi, ...,x n ) be a distribution function 
of N particle. If we want to point out that f(x\, xn) depends on t we will 
write /(xi, xjv|t). Let 

fi(xi) = J dx 2 ...dx N f(x 1 ,...,x n ), (147) 



and 



f 2 (x 1 , x 2 ) = J dx 3 , dx N f(x u x N ) 
be marginal distribution functions. Put by definition 



(148) 



pi(xi) = iV/i(xi), 
p 2 (xi,x 2 ) = ^ 2 /2(xi,^). (149) 

If A is a function on the phase space T, T = M 67V and 

N 

1=1 

then 

(A) = / /(arx, ...,x n )A(xi, ...,x N ) = 



= J dxA{x) Pl (x). (151) 



Now if A is a function on the phase space 

A = ^2A(xi,Xj) (152) 

in the limit of large we find 

(A) = J dx 1 dx 2 p 2 (x 1 ,x 2 )A(x 1 ,x 2 ). (153) 
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Let us introduce also three-particle distribution function: 

f 3 (x 1 ,x 2 ,x 3 ) = J dx 4 ...dx N f(x 1 ,...,x N ). (154) 

Let us derive the equation for f(x). At first let us write equation of motion 
for f(xi, x N ). We have 

o n 

— f(x u x N | t) + ) —V i f{x 1) x n \t) 

at ^ m 

i=i 

dV(qj - qj) df(x 1 , ...,x n \t) 

in 

This equation is only an infinitesimal form of the Liouville theorem. Let us 
multiple this equation by iV and integrate over dx 2 , dx^. Suppose that 
f(xi, xn) is a function of rapid decay of momenta. This assumption admit 
us integrate over pi by parts. We find: 

-^,P\\X\\t) + —Vpiix^t) 
ot m 



+ / dx 2 — ^-p 2 (x u x 2 \t) 
J mdq 2 



^dV^-q^dp^Mt) (156) 
oqi ox i 

Note that we kept here boundary term. Let us now talk about derivation of 
kinetic equation. According to a standard prescription we put pz{x\, x 2 ^x 3 ) = 
in equation for p 2 (xi, x 2 ). We find the following equation for p 2 : 

j t p 2 (x 1 (t),x 2 (t)\t) = 0, (157) 

where (xi(t),x 2 (t)) is a solution of corresponding two-body problem. 

Condition of correlation breaking. In purpose of simplicity we con- 
sider only translation-invariant matter. Usual correlation-breaking condition 
has the form 

p 2 (x u x 2 |0) = h{p' 1 (x 1 ,X2))h{p' 2 (x 1 ,x 2 )). (158) 

Here h is a function on momenta-space of one particle. We consider only 
translation-invariant gas, so h depends only of momentum. 
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Pi and p 2 are momenta of particles 1 and 2 at t = — oo if at t = their 
coordinates and momenta was x\ and x 2 respectively. 
Proposition. 

d 

—p 2 (xi,x 2 ) = 0. (159) 

Indeed, according to (11571) 

p 2 (x 1 ,x 2 \t) = p 2 {x°, xl\o), (160) 

where x\ and x® are phase coordinates of particles 1 and 2 respectively at a 
moment t = 0. Therefore 

p 2 (xi,x 2 \t) = hip'^xl^hip'^xlxl)). (161) 

But the points x® and x® come to the points Xi and x 2 after the time t. So 

biM,^).^^?.^)) = (Pi{xi,X2),P 2 (xi,x 2 )) 7 and 

h(p' 1 (x u x 2 ))h(p 2 (x 1 ,x 2 )) = p 2 (x!,x 2 \0). (162) 

In result 

p 2 (x 1 ,x 2 \t) = p 2 (x!,x 2 \0). (163) 

The proposition is proved. 

It follows from equations (11 5T[) and (11591) that: 

(-V 1 + ^V 2 )/(x 1 ,x 2 |t) 
m m 

,dV(q 1 -q 2 ) d 3V(q l -q 2 ) d 

= o tt- + o ~— / zi, z 2 (164) 

c% 9pi dq 2 dp 2 

The function ^i(p) can be found from the following equation 

pi{x)= lim — I p 2 (x 1 ,x 2 )dx 2 . (165) 

But in zero order of gas parameter the particles are free and 

Pi(x) = h(x). (166) 
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Formula (11641) is usually used for transformation of r.h.s. of equation (11561) 
to the collision integral. From other hand the equation 

d 

—p 2 (x 1 ,x 2 ) = (167) 
at 

shows that there no irreversible evolution in the system. From other point 
of view we will show that the last term in the left hand side of (11561) is equal 
to the scattering integral. 

For simplicity we will show the case v% — 0, v — — . The general case 
can be reduced to this case by means of Galilei transformation. So let us 
consider the integral 

/= lim / d 3 p 2 f dg 2 ^A p ( ,0,p 2 ,g 2 ), (168) 
R->oo J J mdq 2 
v R 

where Vr is a ball of radius R with the center at zero. Let us integrate over 
dq 2 by using Gauss theorem. We find 

/= lim d 3 p 2 dS—cosipp(0,0,p 2 ,q 2 ). (169) 
i?,^oo J J m 

Sr 

Here Sr is a boundary of Vr and ijj is an angle between to rays: first of them 
is parallel to p 2 , second starts from zero and comes throw q 2 . We have 

lim / d 3 p 2 / dS— cos ip x 



R->oo J J m 

Sr 

MPi(0,0))% 2 (p 2 ,g 2 )). (170) 

Let us suppose that the particle scatters only then they are not too far from 
to each other. Then 

h(jp' x {Q,Q))h{p' 2 (jp 2 ,q 2 )) = pi(p 2 )pi(0) (171) 

for all q 2 G Sr\0, where O is a small neighborhood of the point go := ^J-^ 
Sr. Diameter of O is approximately equal to diameter of suppV. Therefore 
the integral / is not equal to zero and equal to 

/ = / d 3 p 2 — [ 2%bdb 
J m J 

x {pM((P2, q 2 (b)), (0, 0)) Pl (p' 2 ((p 2 , q 2 (b)), (0, 0)) 

-pi(P2)pi(0)}, (172) 
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where b := q2 — Qo- But the right hand side of (11721) is a usual collision 
integral. 

Therefore if we keep boundary terms in BBGKI-chain we obtain the ki- 
netic equations without collision integral. 

15 Conclusion 

In this paper we have studied the problem of divergences in Keldysh diagram 
technique which arise if the matter is non-equilibrium. We have considered 
some wide class of divergent diagrams and have proposed a method for renor- 
malization of this divergences. We use this method for renormalization of 
one- and two-chain diagrams. 

A general thesis that we want to illustrate in this series of papers consists 
in follows: to prove that the system tends to thermal equilibrium one should 
to take into account its behavior on its boundary. In the last section we 
have shown that some boundary terms in BBGKI-chain which are usually 
neglected in Bogoliubov derivation of kinetic equation compensate scattering 
integral in kinetic equation. 
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